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A spin transport model is employed to study the effects of spin dephasing induced by diffusion-driven transit- 
time uncertainty through semiconductor spintronic devices where drift is the dominant transport mechanism. 
It is found that in the ohmic regime, dephasing is independent of transit length, and determined primarily by 
voltage drop across the spin transport region. The effects of voltage and temperature predicted by the model are 
compared to experimental results from a 350-micron-thick silicon spin-transport device using derived mathe- 
matical expressions of spin dephasing. 



I. INTRODUCTION 

Employing electron spin in semiconductor devices could 
potentially overcome scaling issues with modern charge- 
based electronics, providing more power efficiency, higher 
speed, and greater functionality ill j&l |3l |4l, |5l y, |7tl . In many 
spintronic device designs, spin precession during transport 
from injector to detector plays a primary role in determining 
device output characteristics. This precession is caused by a 
torque exerted either by efFective[8J or reall,9i, ,10] magnetic 
fields oriented perpendicular to the spin direction. 

Because of magnetic field inhomogeneities or transit time 
uncertainty, transported spins will aiTive at the spin detector 
with a distribution of precession angles. When the angular 
distribution width approaches 2n, the contributions of spin- 
polarized electrons to the detector output is reduced by partial 
signal cancellation of oppositely-oriented spins, in a process 
called "dephasing", or "Hanle effect" lilTl[T2l[T3ll . 

Dephasing in spin-transport devices (especially when due 
to transit-time uncertainty caused by spatial diffusion or path- 
length variation) cannot be reversed, as it can be with spin- 
echo techniques in electron spin resonance (where transport 
plays no role) 1 14]. Therefore, spin dephasing is a physical 
process of fundamental importance to semiconductor spin- 
tronic devices. Here, we present analysis of a spin drift- 
diffusion model to quantify the effects of dephasing in a 
regime where charge transport is dominated by drift, as in 
recent experiments liisl fl6l [TtIi . We examine the apparent 
independence of relative dephasing (precession fringe vis- 
ibility in magnetic field spectroscopy) to injector-detector 
distance (transit length) and quantify the effects of voltage 
bias across the transport region. These conclusions are bol- 
stered through comparison to experiment, using silicon spin- 
transport devices|16|. In addition, the effects of temperature 
on spin dephasing are discussed. 



11. MODEL 

Since spin detection is invariably a projection of the elec- 
tron spin on a fixed measurement axis, the change in output of 
a spintronics device employing precession will be determined 



by cos 6 (to first order), where 6* = w ■ t is spin angle, and oj is 
spin precession frequency and t is transit time from injector to 
detector over transit length L. Clearly, transit-time uncertainty 
At gives rise to precession angle uncertainty lS.9 - oj - At. In- 
creasing precession frequency (w = gfigB/fi, where g is elec- 
tron spin g-factor, //^ is the Bohr magneton, fi is the reduced 
Planck's constant, and B is perpendicular magnetic field) by 
increasing the magnetic field therefore increases the angular 
uncertainty AO and causes dephasing. A complete model for 
dephasing during transport from injector to detector is there- 
fore needed. 

Electron transport in a semiconductor is well-described 
by the drift-diffusion equation. To model spin transport 
accurately, finite spin-lifetime must be included within the 
relaxation-time approximation: 

ds d^s ds s 

where s is the spin density, T^f represents the spin relaxation 
time, D is the diffusion coefficient for electrons, and v is the 
electron drift velocity. In the "ohmic" regime, v - y.E, where 
/i is mobility and E is electric field. Eq. [1] governs the evolu- 
tion of spatial distributions of electron spin, and can be solved 
easily for initial conditions (at f = 0) of an ensemble of spins 
which are all spin-polarized in the same direction at the injec- 
tor (at X - 0). In other words, we find the Green's function 
describing the evolution of a Dirac delta: 

s(x,t) = — -^e "'f . (2) 

The most salient features of this gaussian solution for our 
subsequent analysis are that the center of the distribution trav- 
els with velocity v toward the detector, and that the spatial 
distribution width increases in time as VOf. 

Holding X - L, this solution describes the distribution of ar- 
rival times t at the detector. Therefore, the precession-induced 
change in total device output signal (Ale) comprised from all 
electrons arriving with different transit times < f < oo, each 
contributing cos 6 due to spin precession, is 
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FIG. 1 : (a) Comparison of the spin precession signal output accord- 
ing to Eq. [3] for devices with L = lOO^m (black) and lO/im (light 
grey) transport layer thickness, showing identical dephasing. The 
voltage bias is 2V, mobility is 1400 cm^/V s, diffusion coefficient is 
36 cm^/s, and spin lifetime is 100ns for both simulations. Several 
oscillation extrema are labeled by the average precession angle n, 
2n, 3n, 4n. (b) The independence of dephasing to transport length 
regardless of oscillation order is shown using Eq. [17] calculated for 
lOfi m< L < lOOyU m. 



This integral expression can be easily evaluated to compare 
to empirical observations. A typical experiment consists of 
a spectroscopy where perpendicular magnetic field strength 
B is varied; spin precession oscillations (where extrema cor- 
respond to average precession angles in integer multiples of 
tt) are suppressed in larger magnetic fields due to a larger 
precession-angle uncertainty. This suppression due to dephas- 
ing, or "Hanle effect", can be seen in results evaluated from 
Eq. [3] shown in Fig. More relative dephasing increases 
spin decoherence at the detector and reduces the number of 
oscillations seen. 



ni. DEPHASING ANALYSIS 

Because the integrand in Eq. [3] is dominated by the gaus- 

sian term e , the exponentially decaying part e ^'f can be 
ignored if the spin lifetime r ,/ is suitably long. Then we have 



Mr 



{L-vD- 

-e *Di cos utdt. 



n 1 

Jo 2^^T, 

Making the variable substitution = cot converts Eq. |4]to 
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Application of the following transformation 
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(where A is a transit-length scaling factor) converts Eq. |5]to 
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Other than an overall multiplicative factor of A, Eq. |7]is 
identical to Eq. |5j the spin precession model is therefore 
virtually invariant to this transformation. In terms of device 
physics, ohmic transport (v - iiVIL, where V is the volt- 
age bias on the transport layer) automatically accounts for 
decreasing the drift velocity by a factor of A if the length L 
increases by the same proportion at constant voltage; the first 
two elements of the transformation are therefore satisfied. The 
third element of the transformation in Eq. |6]is accounted for 
by applying a quadratically weaker magnetic field. The ex- 
perimental outcome of this result is that in measurements of 
different devices, we can expect the same number of preces- 
sion oscillations despite variation in transport lengths, assum- 
ing the applied voltage V is identical. 

Again, the preceding assumes that the exponential spin- 
decay term is negligible. This is true when the spin lifetime is 
much longer than the time it takes the gaussian spin distribu- 
tion of width d to enter into the detector (transit time uncer- 
tainty. At = i}. Since for drift-dominated operation 



d^ ^fl>r^ JD- 



(8) 



we have 
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If Tsf » Ar, then finite spin-lifetime can be ignored and 
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FIG. 2: Dephasing (defined by Eq. 117b for the precession oscillation 
extrema calculated with Eq. |3]as a function of applied voltage across 
a 350 fim transport region with mobility /j = 1400cm^/V s (symbols). 
Lines are analytically given by Eq. [T3]with corresponding values of 
6 indicated by the legend. Note exceptional agreement between these 
two complementary means of determining dephasing. 



variations in L make negligible difference to the evaluation of 
Eq. |3] For characteristic values D - lOOcm^/s, v = 10*cm/s, 
and Tsf - 100ns, this condition is satisfied for L <s: 100cm, 
which is certain to hold for semiconductor devices that are 
typically many orders of magnitude smaller. 

Fig. [Ha) shows simulation results for different device trans- 
port layer thicknesses L. The black and light grey curve cor- 
respond to Li = lOO/^m and L2 = AL[ = lOfi m, respectively. 
The simulation is performed using fj. = 1400cm^/Vs, D = 36 
cm^/s and voltage bias V - 2W. The spin lifetime in both sim- 
ulations is 100ns. For the sake of comparison, the y axes are 
slightly shifted relative to each other. It is very clear that the 
spectroscopy shape and precession oscillation extrema num- 
ber is exactly the same, despite the difference in signal mag- 
nitude, which is roughly equal to A (a factor of 10). 

Another way to show the dephasing invariance to transit 
length is to calculate the relative uncertainty in the distribution 
of final precession angle 



A0 oj-At 
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If transport is dominated by drift in the applied electric 
fieldiisj, the transit time is given by 
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(Note the quadratic dependence of transit time on transport 
length at fixed voltage.) Applying ohmic transport v = fiV/L 
to Eq. IHgives 



d = L ^JD/(fiV). 
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FIG. 3: (a) Side-view and (b) schematic conduction band diagram 
of our 350-micron-thick undoped single-crystal silicon spin transport 
device used to compare experimental results to the dephasing model 
presented. 



Since the width of the transit time distribution At is d/v, 
the uncertainty in the distribution of final precession angle 
at the detector is given by application of Eq. [TOl 
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This result (valid at precession extrema when the average 
spin direction is either parallel or antiparallel to the measure- 
ment axis) is independent of the transit length L, so that we 
can expect the same amount of dephasing at the same volt- 
age bias through the transport layer, regardless of the distance 
from injector to detector for any fixed precession angle (as- 
suming ohmic behavior, v - jjE, where E is internal electric 
field). Although Eq. [T3]is independent of transit length L, it is 
clearly dependent on voltage bias V. Specifically, it predicts 
that dephasing is suppressed with an increase in V, which is 
intuitively expected since diffusion will play a smaller role as 
drift becomes stronger in a larger electric field. 

To further study the dephasing effect, it is necessary to es- 
tablish a proper quantitative empirical standard for compar- 
ison to experimental measurements not dependent on data 
fitting using evaluation of Eq. |3] If drift dominates and 
L/ V » At, the factor proportional to 1 / V? in Eq. |3] can 
be approximated as a constant. Then, (assuming again that 
Tsf » At) the distribution of precession angles is gaussian: 



(14) 



so we have 
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We can use this analytic result to find the angular distribu- 
tion width from Ale spectroscopy measured empirically. The 
ratio of the magnitude of the central maximum for B - Q (cor- 
responding to zero precession: 6-0 and hence A9 = 0) to the 
signal value for extrema at 6 = nn (n = 1, 2, 3...) is therefore 
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It is therefore possible to compare the dephasing effect in 
different experiments by studying the signal magnitude de- 
crease of the same order extrema, without having to fit to the 
complete model (Eq. O. Fig. [TJb) shows the distribution 
width AO for the several precession oscillation extrema based 
on Eq. [IT] for simulations of devices with different transport 
layer thickness from lOyum to 100/im. The same voltage bias, 
2V, is applied on the transport layer It is obvious that, de- 
spite an order-of-magnitude change in the transit length L, the 
dephasing is virtually constant, further confirming the conclu- 
sions drawn from Fig. ^a). 

Voltage-controlled dephasing in a device with L = 350/im 
from our simulation is shown using Eq. [T7] in Fig. |2] 
(symbols). At fixed voltage, the changes in dephasing for 
successive extrema are approximately equal, because AO - 
g/UBBAt/h and the magnetic field period of precession oscil- 
lations is approximately constant. We compare the dephasing 
behavior of several precession oscillation extrema from these 
calculations to the behavior predicted by Eq. [13] (solid lines). 
Despite there being no free parameters for fitting, we find ex- 
cellent graphical agreement between Eqs. [13] and [IT] for all 
precession oscillation extrema examined, confirming the va- 
lidity of our obtained expressions. 
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FIG. 4: (a)-(e): Experimentally-measured spin precession spec- 
troscopy at different accelerating voltages lOV, 14V, 16V, 18V, and 
80V, respectively, using spin transport through 350 micron-thick un- 
doped single-crystal Si at 150K. Unprecedented spin coherence is 
evident in (e), showing at least 15 full spin rotations. 



IV. EXPERIMENTS 

We now compare these results to experimental data, us- 
ing all-electrical hot-electron spin injection and detection 
to study spin transport through 350-micron-thick undoped 
single-crystal Si(100)|[T3l- The device scheme and corre- 
sponding band structure are shown in Fig. [3 a) and (b). Trans- 
port proceeds from top to bottom: equilibrium spin-polarized 
electrons from the Co84Fei6 cathode tunnel through the AI2O3 
barrier, becoming hot electrons in the nonmagnetic Al/Cu an- 
ode thin film. The electrons coupling with conduction band 
states over a Schottky barrier in the silicon transport layer on 
the other side forms the injected spin polarized current. Drift- 
dominated transport in the undoped silicon transport layer is 
controlled by Vc\, the "accelerating voltage". After trans- 
port through this silicon layer, the electrons are ejected from 
the conduction band and into the second ferromagnetic metal 
(Ni8oFe2o) and the ballistic component of this current is col- 
lected by an n-type silicon substrate below, forming the spin- 
transport signal Ic2- This current is dependent on the pro- 
jection of electron spin direction onto detector magnetization 
axis of the Ni8oFe2o layer|16]. A magnetic field aligned per- 



pendicular to the spin direction and parallel to internal elec- 
tric field (and hence drift velocity) induces spin precession 
during transport from injector to detector. (This device is es- 
sentially the solid-state - and Si-based - version of the first 
bare-electron g-factor experiment which used spin precession 
during ballistic transport of high-energy electrons in vacuum 
and Mott scattering for spin filteringi 19].) 

Fig. |4|a)-(e) shows precession measurements at 150K using 
this device with accelerating voltage lOV, 14V, 16V, 18V and 
80V, respectively. As the accelerating voltage increases, more 
oscillations are visible; at the highest voltage applied, more 
than 15 full spin rotations are evident by the number of oscil- 
lation extrema in Fig. |4|e). The length independence derived 
earlier (Eq. [T3] ) explains the high degree of spin coherence in 
these precession measurements, despite the 350 microns of Si 
separating injector from detector 

Using Eq. [T7J we plot the dephasing at 150K as a function 
of voltage for several precession oscillation extrema in Fig. 
[5] and compare to the y '^^ behavior predicted by Eq. [13] 
Although the empirical behavior is qualitatively correct, the 
remaining discrepancy is likely due to our model assumption 
of ohmic transport, which at these voltages is not precisely 
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FIG. 5: Dephasing (calculated from experimental data with Eq. I17t 
for several precession oscillation extrema with different voltage bias 
across the 350-micron-thick undoped single-crystal Si transport layer 
at 150K. The model-predicted V"''^ behavior is also plotted for com- 
parison (dotted line). 
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FIG. 6: (a)-(c) Temperature dependence on spin dephasing for our 
Si spin transport devices at three temperatures 20K, 60K, and lOOK, 
respectively. Panel (d) plots the dephasing (calculated using Eq. I17t 
for 20K<T<100K, showing the unexpected non-monotonic depen- 
dence unaccounted for by application of the Einstein relation to our 
model. 

upheld due to onset of velocity saturation ifisl . 

The Einstein relation (or "fluctuation-dissipation theo- 
rem"), valid for non-degenerate charge density, dictates 
D/jj - ksT/q, where ks is Boltzmann's constant and T is tem- 
perature. Because our injected cuiTents are low (^ 1 - lOyuA), 
and the spin-injector area is large 50 x 100/im^), our ex- 



perimental electron density n = ss lO'^ cm^ is clearly 
non-degenerate in Si which has an eff'ective conduction DOS 
» 10'''cm''. If the Einstein relation is also valid for the non- 
equilibrium electrons we generate using ballistic electron in- 
jection with our tunnel-junction emitter, Eq. [T3] implies that 
the relative dephasing AO/ 6 oc Vt. Therefore, more dephas- 
ing should be evident at higher temperature. 

Contrary to this expectation, we find a non-monotonic be- 
havior of experimental dephasing as a function of temperature 
for experimental measurements at fixed accelerating voltage, 
as shown in Fig. |6ta)-(c). The dephasing is small at the lowest 
temperature (20K), then rises to a maximum at approximately 
60K, and drops gradually toward lOOK, as shown in Fig. I^d). 
Clearly this observation is not compatible with the Einstein re- 
lation, and a more sophisticated model than the one presented 
here is required to capture the correct behavior It is interest- 
ing to note, however, that the sharp decrease in dephasing at 
low temperatures corresponds coincidentally to the onset of a 
negative differential mobility (NDM) regime in Si(lOO) below 
~40K (Ref. |20]), which in other semiconductors gives rise to 
the Gunn effect: transport is via soliton-like spatial charge do- 
mains that suppress longitudinal diff'usion and would therefore 
constrain spin dephasing in these spintronic devices. How- 
ever, NDM occurs in Si(lOO) only in a window of electric field 
values (a!50-150V/cm) that are substantially smaller than the 
ones used here. It is more likely that the observed reduction of 
dephasing is the result of the consequences of (unintentional) 
doping impurity freeze-out. 



V. CONCLUSION 

Using a spin precession model based on the Green's func- 
tion of the spin drift-diffusion equation, we determined the 
expected effect of transit length and voltage changes on spin 
dephasing. Most importantly, we determined that for long 
spin lifetimes, the relative spin dephasing is independent of 
the transit length. In future spintronic devices, this length in- 
dependency may enable multiple gate operations and long dis- 
tance transport in precession-dependent devices as long as the 
finite spin lifetime allows. 

We also constructed an empirical measure of spin dephas- 
ing and used it to characterize both the results of simulation 
and experiment. Dephasing is shown to be inversely depen- 
dent on the square-root of voltage drop that drives drift, both 
theoretically and experimentally. 

For model results closer to experiment, the true 3- 
dimensional device geometry, magnetic fringe fields from in- 
jector and detector contacts, and non-ohmic transport, (all of 
which will give other systematic sources of dephasing) should 
be taken into account. 

This work was supported by DARPA/MTO and ONR. 
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